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Abstract 



Shapes and high-spin properties of nuclei from the neutron-rich (A^>56) 
zirconium region are calculated using the Nilsson-Strutinsky method with 
the cranked Woods-Saxon average potential and monopole pairing resid- 
ual interaction. The shape coexistence effects and the competition between 
rotationally-aligned I/111/2 neutron and 1^9/2 proton bands is discussed. Pre- 
dictions are made for the low-lying superdeformed bands in this mass region, 
characterized by the intruder states originating from the J\f=5 and 6 oscillator 
shells. 

PACS number(s): 21.10.Ky, 21.10.Re, 21.60.Jz, 27.60.+j 
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I. INTRODUCTION 



Nuclei from the heavy-Zr region (Z~40, A^>56) exhibit a wealth of structure phenomena 
such as, shape coexistence effects, dramatic variations in quadrupole collectivity, strong oc- 
tupole correlations, and the existence of low- lying intruder states Most nuclei from this 

region are on the neutron-rich side and, therefore, they are not accessible using the standard 
in-beam techniques. The heaviest known isotopes in the A~100 mass region, namely ^''^Y, 
^°''Zr, ^^''Nb, ^^^Mo, ^^^Tc, and ^^^Ru, have recently been observed ||^ using the technique of 
in-flight isotopic separation of projectile fission fragments at relativistic velocities. Proper- 
ties of /3-decay of ^^^y, ^^^Zi^ ii^Nb, i^°Mo, ^^^Tc, and ^^^Ru, have been studied §,§ at the 
IGISOL facility. These nuclei determine the present experimental spectroscopic boundary 
of this region. 

A powerful tool to carry out spectroscopic studies in the neutron-rich systems is to 
analyze the prompt gamma rays from nascent fission fragments. In such measurements it 
is possible to approach spectroscopically very neutron-rich nuclei such as ^'^^Sr, ^°^Zr, ^°^Mo 
|0J^, and ^^^Ru 0. Rotational structures in this region have also been studied using breakup 
fusion reactions |T^, (t,p) reactions |Tl|, heavy-ion-induced fusion []12|, incomplete fusion 
T3| , neutron- induced fission [ll4| , p!5| , and deep inelastic reactions |[l6|jr^ . 

In spite of many experimental efforts, detailed information on high-spin properties of 
y4~100 nuclei is still scarce. The first systematic study of high-spin states in even-even 
nuclei, as well as odd-A nuclei in this region, was reported in Ref. In the even-even Zr 

in-ii 



and Mo isotopes, the yrast lines are known up to I'^ 

108 



--12- 



In several nuclei such 



as io2Ru 13, j^u p 



16| , a crossing between the ground band 



"2Ru H, and Mo 

g-band) and the aligned band (s-band) has been observed. Also, in a number of Ru and Mo 
side-bands have been studied 
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Zr m, lo^zr m 



T|,|2T]-||]. 
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High-spin states in odd-A nuclei have 
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Nb [gH, i°3Mo [gg, ^o^Mo |0 



32| . For the heavy Sr isotopes, spectroscopic studies ||33| , |34 



are 



isotopes 

been found in ^°^Sr 
i03Rh p, and 109'iiiRu 
still in a preliminary stage, and practically nothing is known on the neighboring Kr and Se 
nuclei. It is conceivable that with the help of the new generation gamma-ray arrays such as 
GammaSphere and EuroBall, and with the help of radioactive beam facilities, the borders 
of the heavy-Zr regions will soon be expanded towards even heavier systems. 

The strong dependence of observed spectroscopic properties on the number of protons 
and neutrons, makes the neutron-rich A^lOO nuclei a very good region for testing various 
theoretical models. According to the mean-field-based calculations, strong shape variations 
in this region may be attributed to shell effects associated with large spherical and deformed 
subshell closures in the single-particle spectrum. The single-particle diagram, representative 
of the nuclei from the heavy-Zr region, is shown in Fig. |I]. The strongest shell effects are 
expected for a spherical shape with the magic numbers 28, 50, and 82, and also for the 
subshell closures at Z=40 and A^=56. The transition to a deformed region occurs near 
A^=58 and 76; the strongest shell polarization towards deformed shapes is predicted for 
A^=60-72. Finally, the best candidates for superdeformed shapes are nuclei with particle 
numbers near 42 and 64. [Shell-correction landscapes for the Zr- region are shown, e.g., in 
Ref. m 
model).] 



(folded- Yukawa model), Ref. |]36| (Woods-Saxon model), and Refs. ||37| , |38[| (Nilsson 
According to calculations based on the mean-field approach, the occupation of I/111/2 
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neutron and lgg/2 proton orbitals is essential for understanding the deformed configurations 
around ^°°Zr |^9|. The best examples of shape coexistence in this region are Sr, Zr, and Mo 
isotopes with iV~58. In the language of the deformed shell model, the onset of deformation 
around A^=58 can be associated with the competition between the spherical gaps at Z=38, 
40, and A^=56, and the deformed subshell closures at /52~0.35 at particle numbers Z=38, 
40, and A^=60, 62, and 64. Theoretically, the delicate energy balance between spherical and 
deformed configurations depends crucially on the size of these gaps in the single-particle 
energies. As discussed in Refs. |pO|-^ the deformation onset at N^58 results from the 
subtle interplay between the deformation-driving neutron-proton energy (varying smoothly 
with the shell filling) and the symmetry-restoring monopole energy responsible for the shell 
effects. 

Equilibrium deformations and moments, potential energy surfaces, microscopic structure 
of coexisting configurations, and shape transitions in the heavy-Zr region have been calcu- 
lated by many authors In most cases, calculations show large deformations in Sr, 
Zr, and Mo isotopes with A^>60. However, the details of the shape transition near A^=58 is 
predicted differently by various models, the onset and rapidity of this transition being very 
sensitive to the model 0. 

One of the most interesting features of nuclei from the y4~100 region is the richness of 
various structural effects that occur at high angular momenta. Many of these effects have a 
straightforward interpretation in terms of the interplay between deformation, pairing, and 
the Coriolis force. In the deformed nuclei near ^°^Zr, the alignment pattern is expected to 
be rather simple. While the first proton crossing can be associated with the breaking of the 
1(79/2 pair, the lowest neutron crossing is due to the l/iii/2 alignment. 

For even-even systems, band crossings have been seen in the g-bands of ^°^^^^^Ru the 
A^=60 isotones ^^^Ru ||JTg,|TB|JT3] and ^^^Fd the A^=58 isotones i°°Mo [|,|Tg], ^°^Ru Jig, 



and Pd p6[, and the A^=56 isotone Ru [^. For proton numbers Z>46 the interplay 



between lgg/2 protons and I/111/2 neutrons is well established [^,^]. In the case of ^°^Ru, 
the 1/111/2 neutron crossing is supported by blocking arguments in neighboring odd-A^ and 
odd-Z nuclei |10,31]. Namely, in the A^=58 nucleus ^°^Rh there is a clear crossing in both 
7r=— and 7r=+ 1-q.p. proton bands, as well as in the 1-q.p. 7r=+ bands in ^''^'^^^Ru. 
However, no crossing has been observed in the vr=— , l/iii/2 bands in ^''^'^^^Ru. Similarly, 
the presence of a band crossing in the yrast line of even-even nuclei ^°^Ru and ^°^Mo ||16|| , 
and its absence in the 7r=— 1-q.p. bands in the neighboring odd-A^ systems, suggest that 
l/iii/2 neutrons are involved. 

In addition to blocking arguments, there exists other experimental evidence for the pres- 
ence of the l/iii/2 neutron unique-parity states in this mass region, notably the direct obser- 
vation of negative-parity, very collective bands in ^^^'^^^Zr, assigned to the [532 5/2] (1/^11/2) 
Nilsson orbital [0,11111 • 

In the Ru and Mo isotopes with A^>60, the role of (static or dynamic) triaxial shape 
degrees of freedom has been discussed. In ^''^~^°^Mo and ^''^"^^^Ru, the excitation energies of 
the 2^ bandheads are low. Furthermore, in the ^"'^'^^^Mo isotopes, the band built upon the 
excited /'^=4"'" state, a candidate for the double-7 vibrational excitation, has been observed 
p^ , |23[] . On the other hand, based on the energy differences and the E2 branching ratios 
analysis, triaxial deformations have been suggested for ^°^~^^^Ru [jl4[. The development of 



triaxiality at high spins has been addressed by measurements of 7-band— >g-band branching 
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ratios M and by lifetime measurements |T9|,pO 



Little is known about shapes of lighter elements in this region. Recent laser measurements 
of isotope shifts in the heavy Kr isotopes 0,0] indicate the absence of a sudden shape 
transition around N=60. 

The main objective of this study is to investigate the importance of shape changes and 
pairing correlations in the description of high-spin bands in nuclei near ^°^Zr. Predictions 
have been made for those nuclei which show the best prospects for triaxial shapes and for 
superdeformed bands at relatively low excitation energies. 

The paper is organized as follows. The Woods-Saxon model and the total routhian 
surface technique are described in Sec. 0. The results of ground-state calculations in the 
y4~100 region are given in Sec. |T|. In particular, ground-state equilibrium deformations 
and their dependence on the treatment of pairing correlations and macroscopic energy are 
discussed. Cranking calculations are presented in Sec. together with the analysis of typ- 
ical alignment patterns, band structures, and shape changes. The summary and conclusions 
are contained in Sec. 0. 



II. THE MODEL 



The macroscopic-microscopic method employed in this work is an approximation to the 
HF approach [p2| , |93| . Its main assumption is that the total energy of a nucleus can be 



decomposed into two parts. 



E = E„ 



(1) 



where -Emacro is the macroscopic energy and -Emicro is the microscopic energy (shell correction) 
calculated from a non-self-consistent average deformed potential. 

The shell corrections were obtained using the deformed Woods-Saxon (WS) potential 
of Ref. [0. For the macroscopic energy, two energy formulas were employed: the liquid- 
drop (LD) model of Ref. and the Yukawa-plus-exponential model (finite-range liquid- 



drop (FRLD) model) of Refs. |]96|j97[| . The latter model is softer to shape distortions and, 
consequently, is expected to favor more deformed nuclear shapes |98 Hl00| • 

The nuclear surface was defined by means of the standard multipole expansion. 



R{n) = c{a)Ro 



(2) 



with c{a) being determined from the volume-conservation condition and Rq = tqA^^^. [Axial 
deformations Px are defined as Px = axo-] 

In our analysis of equilibrium shapes we did not consider the odd multipolarity deforma- 
tions. We are mainly interested in the dominant features related to the leading quadrupole 
effects. Secondly, as discussed previously [p^ , p.01|| , the static octupole deformations are not 
expected to be present in well-deformed neutron-rich 74~100 nuclei. 

The shell correction was calculated according to the prescription given in Ref. |Q . The 
results of the calculations presented in this paper were obtained with a smoothing width of 
49.2/^4^/^ MeV and a sixth-order curvature correction. All single-particle states lying below 
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a cut-off energy of 131.2/^1/3 MeV above the Fermi level were included. Our calculations 
are based on a three-dimensional minimization on the mesh, and the variation of the average 
potential with Z and N has been suppressed (i.e., the calculations were performed using 
single-particle energies of a representative central nucleus). The related errors in energy 
surfaces introduced by this procedure are known to be small [|102| . 



In the /=0 calculations, pairing energies were computed using either the standard BCS 
treatment or the particle number projection method before variation (PNP). The latter 
method is expected to be more precise and reliable in the regions of low single-particle level 
densities where the BCS treatment breaks down. For the pairing calculations we took the 
lowest Z OT N single-particle orbitals for protons or neutrons, respectively, counting from the 
bottom of the potential well. The pairing strengths were determined using the average gap 



method. The optimal choices for the average gap parameterization, discussed in Ref. [ 103 
are adopted for this study. Namely, in the PNP variant we used average gaps of A„ = 
8.6/y/AMeV and Ap = 9.9/y/AMeV, and in the BCS variant we used A„ = 13.3/y/AMeV 
and Ap = 13.9/71 MeV. 



The parameters of the WS potential were taken from Ref. however, they were 

adjusted to describe the rapid shape transition between ^^Zr and ^''^Zr. These parameters 
differ from the standard set ||104|| only in the values of the spin-orbit potential strength. 



Namely, we adopted the values Aso,n=39 and Aso,p=32 instead of the standard values Aso,n=35 
and Aso,p=36, respectively. 

In order to study the variations of potential energy surfaces at /=0 (Sec. PD , three sets 
of calculations were performed: 

VI Axial calculations in a three-dimensional deformation space (/32, /S^, (3^). They were 
carried out by setting up a three-dimensional mesh with the mesh steps equal to 0.05 
0.04, and 0.04 for (32, P'i, and Pq, respectively. In these calculations we employed both 
LD and FRLD models for the macroscopic energy and both pairing models, i.e., BCS 
and PNP. 

V2 The explicit energy minimization in a three-dimensional deformation space {(32, (^i, (3q) 
performed for each isotope separately (here no assumption about a "central nucleus" 



has been made). In these calculations, we used the set of WS parameters of Ref. ^6 
the FRLD macroscopic energy, and the Lipkin-Nogami (LN) approximation to the 
PNP method (with the same pairing strengths as for the PNP mesh calculations). 



The details of the LN treatment are given in Ref. [|102|] 



V3 Triaxial calculations in a three-dimensional deformation grid {(32, 7, (34). The mesh con- 
sisted of 10 points [0.05 (0.05) 0.50] in the /?2Cos(7 + 30°) direction, 12 points [-0.20 
(0.05) 0.35] in the P2 sin(7 + 30°) direction, and 4 points in the direction of P^. To 
optimize the calculations, the /34-grid was defined relative to the value of (3^ minimizing 
the LD energy at given P2 and 7, /94,ld(/325 7) [|105|| . Namely, the hexadecapole mesh 



was defined as /34,ld(/3, 7) + A/34, A/34 being equal to -0.04 (0.04) 0.08. The hexade- 
capole deformation was included in a way which guarantees the modulo 60° invariance 
of the nuclear shape in the {(32, "y) plane, i.e., the nonaxial hexadecapole deformations 



are present for 7 not equal to multiples of 60° P3]. In this variant we employed the 



LD model for the macroscopic energy and both pairing models, i.e., BCS and PNP. 
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The high-spin behavior (/>0) of y4~100 nuclei was investigated using the total routhian 
surface (TRS) approach of Refs. [|106| , |107 ]. The TRS cranking calculations were performed 
on a three-dimensional deformation mesh including /32, 7, and /?4 deformations. The mesh 
was the same as that described in variant V3. The calculations were carried out for all 
nuclei with Z=38-44 and A^=58-66, and for rotational frequencies ^ti;=0.075 (0.075) 1.50 
MeV. The total routhians were defined as: 



uj=0 



H' 



uj=0 



(3) 



where £''^=*^ is the energy calculated at a given deformation according to the shell-correction 
method, H'^=H-ujjx is the total routhian including the pairing term, and \^>'^) is the in- 
dependent quasi-particle vacuum at a given rotational frequency. Due to the size of this 
computation, the LD model was employed for the macroscopic energy. In our analysis of 
rotational bands, we limited ourselves to the vacuum configurations of even-even systems. 
That is, only the positive-parity, even-spin (signature r=l) bands were considered. 

The average BCS pairing gap was parameterized in a way which facilitates separate 
calculation for neutrons and for protons: A„ = 4.93iV~^/^ and Ap = 4.70Z~^/^ ||103||. For 



nonzero spins the pairing gaps were parameterized as functions of rotational frequency: 




(4) 



The parameter hujc was fixed at 0.75 MeV for neutrons and 0.9 MeV for protons. The proton 
and neutron chemical potentials of the RBCS equations were adjusted at each deformation 
and rotational frequency separately so as to satisfy the particle number equation. 



III. GROUND-STATE EQUILIBRIUM DEFORMATIONS 

The potential energy surfaces (PES) for the neutron-rich even-even Kr-Pd isotopes ob- 
tained in the axial variant VI are shown in Figs, ^and ^ The results of triaxial calculations 
(variant V3) are displayed for selected nuclei in Figs. | (^oO'^o^Zr), | (lo^.ioezr)^ | (100,102^0), 

^ (104,106Mo), 1(104,1061^^)^ g^j^^ |(108'"0Ru). 

The calculated equilibrium deformations obtained by means of the explicit minimization 
at axial shape (variant V2) are shown in Table |, and Table displays the calculated 
equilibrium deformation obtained by a minimization in the (/52, 7)-plane (variant V3). 



A. General properties of potential energy surfaces 



The modification of the WS parameters according to Ref. has some impact only 
on the isotopes with A^=56 and A^=58; the optimized parameters tend to lower the PES 
around the spherical point. The most pronounced differences are obtained for the Sr and Zr 
isotopes where the equilibrium deformations are drastically changed. For other iV=56, 58 
isotones the equilibrium deformations are only slightly modified. More pronounced is the 
dependence on the choice of the macroscopic energy and on the pairing treatment. Namely, 
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the most pronounced minima and the largest equihbrium deformations were obtained in the 
FRLD/PNP model (Figs. |^ and |^). The LD/BCS variant yields rather shallow minima, 
especially for the Mo and Ru isotopes, and the smallest equilibrium deformations. 

For most nuclei considered, the calculations predict two minima corresponding to prolate 
and oblate shapes. In the Kr, Sr, and Zr isotopes, these minima are separated by a relatively 
high spherical barrier. The largest prolate deformations cluster at the middle of the shell, 
i.e., around Z=38, 40, and A^=60, 62, and 64. As discussed above, these particle numbers 
can be associated with large prolate gaps in the single-particle spectrum (see Fig. |l|). 

The barrier height (and the oblate-prolate) energy difference decreases with Z. As will be 
discussed below, in many cases the secondary minima are unstable to triaxial deformation. 
In the PNP calculations, there appear local substructures in the potential energy curves; 
the PES's calculated in the BCS variant are smoothed. The largest effects occur in the Mo 
isotopes where the prolate minimum splits into two components, and in the Ru isotopes 
with A^=72-76, where secondary minima with large deformation ~0. 45-0. 50 appear. 

In triaxial calculations (V3), the Sr and Zr isotopes are predicted to be axial, while the 
Mo-Pd isotopes are calculated to be 7-soft. This is illustrated in the potential energy surface 
(PES) described in Figs I, |(Zr), |, 0(Mo), | and |(Ru). In the PNP variant, the PES's 
of the Mo isotopes with A^=58, 60 show very shallow, 7-unstable valleys corresponding to 

ranging from 0.22 to 0.30. The heavier Mo isotopes with A^=62-66 have shallow triaxial 
minima at 7^20° and (32 ~0.32. All Ru isotopes have been found to be triaxial in their 
ground states, with shallow minima at 7;^20°. Again, the minima calculated in the PNP 
variant are more pronounced than those in the BCS variant. (For more discussion on the 
influence of the PNP on the rigidity of the PES minima, see Refs. |61] , |108|| .) 



B. Equilibrium deformations and quadrupole moments 

The calculated equilibrium deformations obtained by means of the explicit minimization 
at axial shape (variant V2) are displayed in Table |, together with the predicted oblate- 
prolate energy differences, and calculated and experimental ||109|| intrinsic quadrupole mo- 



ments. The oblate-prolate energy difference is defined as AEop=Eo-Ep, where corre- 
sponds to the minimum with negative (32 while Ep corresponds to the minimum with zero 
or positive (32- The quadrupole moments were extracted from the calculated equilibrium 
deformations using the procedure outlined in Ref . ||102|| . In Table |n|, the calculated equilib- 



rium deformations obtained by means of the minimization in the {(32, 7) plane are displayed, 
together with the calculated transition quadrupole moments, defined as [|110|] 

Q(7) = ^Qocos(7 + 30°), (5) 

where Qo=(5(7 = 0°). 

It is instructive to compare the equilibrium deformations obtained in our work with 
those of Chasman calculated using the LD model and the shell correction approach with 
WS potential |Q , those obtained by Moller et al. in the shell correction approach with the 



folded- Yukawa single-particle potential and the finite-range droplet-model (FRDM) ||7J , and 
with the results of the extended- Thomas- Fermi model (ETFSI) of Ref. 
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The heavy Se isotopes are predicted to be oblate-deformed with /32~-0.28 and very large 
values of (3^^0.10. The prolate-deformed minima appear 0.4-1.4 MeV higher in energy. This 
result is consistent with the ETFSI calculations which yield oblate shapes with /?2=^0.31 
for all Se isotopes with iV>58. On the other hand, FRDM favors prolate shapes with 
0.29</32<0.33. 

The heavy Kr isotopes are predicted to be oblate-deformed (/?2~-0.31, /56~0.10) up to 
A^~66. For the Kr isotopes with A^>64, AEop is rather small. This is due to the large oblate 
gap with Z=36. This result is consistent with the ETFSI calculations which yield an oblate- 
prolate shape transition at N^68. The FRDM favors well-deformed prolate ground states 
with /52~0.33. The recent laser-spectroscopy measurements 0,0] show that the observed 
(r^) does not drastically change at A^=60, indicating that there is no sudden increase of 
1/32 1 for A^<60 nuclei. Since the deduced P2 is ~0.3 for A^=60, all these calculations are 
consistent with this experiment. 

The heavy Sr isotopes with 60<A^<72 are predicted to have well-deformed prolate ground 
states in all the models discussed. The calculated quadrupole deformations P2 range from 
0.30 to 0.40. Contrary to the FRDM results, both our calculations and those of the ETFSI 
model yield a rapid shape transition at N^56. 

The heavy Zr isotopes with 60< A^<72 are predicted to have well-deformed prolate ground 
states (0.33</32<40) in all the models discussed. At N^Ti, the transition to oblate shapes is 
predicted. This transition appears much earlier, at A^=66, in the calculations by Chasman. 
The rapid shape change between N=56 and A^=60 is reproduced by our model, by Chasman, 
and by the ETFSI model. In the FRDM, the shape change is calculated to be more gradual; 
the nucleus ^^Zr is predicted to be prolate-deformed with /?2=0.2 (see the discussion in 
Ref. |35] regarding the single-particle structure of the folded- Yukawa potential near ^^Zr). 
As seen in Figs. ^, the heavy-Zr isotopes with A^>66 become 7-soft. 

In the heavy Mo isotopes, the calculated oblate-prolate energy difference decreases, hence 
the competition between low-lying coexisting minima is more pronounced. In contrast to 
the Sr and Zr isotopes, the transition from spherical to deformed shapes in the Mo isotopes 
at iVf»58 is predicted to be gradual, in agreement with experimental findings. The transi- 
tional nucleus ^'^^Mo is very soft and the isotopes io4,io6,io8j^q 

are calculated to have triaxial 

ground-state minima. Similar results were predicted by Chasman, in the Nilsson-Strutinsky 
calculations of Ref. ||5^, and in the Skyrme-Hartree-Fock study The prolate-to-oblate 
shape transition is expected at A^;^68-70 in all these models. 

In the heavy Ru isotopes, the prolate-to-oblate shape transition is predicted at A^?»66 
in our axial calculations and in the FRDM, but much earlier, at A^^60, in the ETFSI 
model. Again, as in the Mo case, our calculations for the Ru isotopes predict 7-softness or 
triaxiality. Similar results were obtained by Chasman, by the Nilsson-Strutinsky calculations 
of Ref. 1^ , and by the Skyrme-HF study . 

Stable triaxial shapes in ^'^®~^^°Ru have recently been reported in Ref. Our calcu- 
lations are consistent with the experimentally deduced 7 values, namely 23° for ^°^Ru and 
24° for ^^°Ru, but fail in reproducing the experimental quadrupole moments. That is, after 
correcting for the effect of triaxiality using Eq. (j^), theoretical quadrupole moments for Mo 
and Ru isotopes are predicted to be too large. Interestingly, axial calculations reproduce 
measured values rather well (cf. Tables | and HD, but they overestimate the experimental 
quadrupole moments for other systems predicted to have near-axial shapes. A possible source 
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of this discrepancy is the shape mixing phenomenon which gives rise to the fragmentation 
of the experimental low- lying B{E2) strength. 



IV. HIGH-SPIN CALCULATIONS 
A. Total Routhian Surfaces: Examples 



Figures |T0|, |Tl], and |T2| show typical examples of total Routhian surfaces for nuclei from 
the heavy-Zr region. The nucleus ^°°Sr, Fig. |10|, represents a well-deformed collective rotor. 
The prolate minimum, /?2~0.35, does not change significantly with angular momentum up to 
J~36. Only at a very high rotational frequency does this minimum become triaxial (7>0°) 
due to alignment of several high-j quasi-particles (see Sec. [iV Bj ). The coexisting oblate 
minimum (7^-60°) is also fairly stable with I. 

In ^°^Zr (Fig. |ll]) the situation is different. There is a shape change from the well- 
deformed prolate minimum to a triaxial minimum with 7<0° associated with the first band 
crossing. The local minimum at 7~-60°, seen in the ^ti;=0.3 MeV surface, quickly be- 
comes non-yrast with increasing uj. At very high rotational frequencies, the superdeformed 
(/32~0,4), slightly triaxial, minimum becomes yrast at 1^50. 

The nucleus ^°^Ru (Fig. O) is representative of a 7-soft system. Here, the competition 



between several near-yrast triaxial minima produces dramatic shape changes. Although 
the triaxial minima are rather shallow at u=0 (see Fig. |^), they become stabilized at high 
spins due to quasi-particle alignment. The superdeformed minimum with /?2~0.5 appears at 
exceptionally low spins. At huj=0.6MeV, this minimum lies only ^2 MeV above the triaxial 
yrast line. 



B. Shape Changes at High Spin 

The systematics of calculated equilibrium /?2 and 7 deformations for even-even nuclei 
with iV=58-66 are displayed in Figs. [T| (Sr), |T| (Zr), |l| (Mo), and |16] (Ru). All the local 
minima that appear in calculated TRS's are visualized as /-dependent trajectories in the 
(/?2,7) space. In the following discussion, rotational bands and their underlying intrinsic 
configurations are grouped into several categories, according to their deformations and in- 
truder content. In most cases, the dependence of equilibrium deformations on configuration 
and angular momentum can be simply explained in terms of two factors: (i) the deformation 
softness of the PES at /=0; and (ii) the position of the Fermi level. A, relative to an intruder 
high-j shell. In this context, recall that the routhian of the strongest-aligned quasi-particle 
has a minimum as a function of 7 that is strongly dependent on the position of A within 
the shell. Namely, if the Fermi level lies at the bottom of the shell, positive values of 7 are 
favored, for a half-filled shell 7~-30° is preferred, and for A in the upper half of the shell. 



the aligned high-j particle drives the core towards 7<-60° Of course, the total 

shape polarization exerted by aligned quasi-particles results from an interplay between the 
contributions of protons and neutrons. 
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1. Prolate shapes with P2 ~ 0.27-0.35. 

The well-deformed prolate shapes are typical of ground-state rotational bands in Sr and 
Zr nuclei. As seen in Table |, ground-state deformations gradually increase from /52~0.28 in 
A^=58 isotones to /52~ 0.33-35 in the heavier systems. 

For all Sr isotopes except ^°^Sr (A^=66), the deformation of the aligned (1^/111/2)^ 2-q.p. 
configuration is nearly the same as that of the ground state band. Indeed, the calculated 
values of P2 are only slightly reduced and the values of I7I are close to zero. For ^°^Sr, the 
neutron Fermi level is higher in the l/iii/2 shell, i.e., between [532 5/2] and [523 7/2] orbitals 
(Fig. 1^). Consequently the (z//iii/2) alignment is reduced and the z//iii/2 crossing is predicted 
at a rotational frequency close to that of the {ngg/2) crossing. 

The values of 7 calculated for the (z^/iii/2)^(7r(y'9/2)^ 4-q.p. configurations depend on 
(see Fig. |T3|). This can be explained in terms of the polarizing effect of the aligned high-j 
particles. Indeed, for N=58 (bottom of the i'hii/2 shell), 7 deformation of a 4-q.p. band is 
calculated to be positive, and with increasing N the system is driven to negative 7-values. 
Interestingly, this polarizing effect is not observed for the (z^/iii/2)^ bands. This can be 
explained in terms of a rather 7-rigid potential of the vacuum configuration (g-band). The 
7-driving force exerted by the aligned (1^/^11/2)^ pair is not sufficient to polarize the system 
towards 77^0. However, the additional deformation-driving effect of the aligned 1(79/2 protons 
is sufficient to make the system triaxial. As expected, the alignment of the second pair of 
l/iii/2 neutrons gives rise to a reduction of P2 in the ('^/iii/2)'^(7r5'9/2)^ 6-q.p. configuration. 
The exception to this pattern is ^^Sr where the increase in j32 can be attributed to the 2/(77/2 
alignment. 

The Zr isotopes are 7-softer than the Sr isotopes (cf. Figs. |10| and |llj), and shape 
changes at high spins are more pronounced. The g-bands are predicted to be prolate. The 
deformation change is associated with the I/111/2 neutron alignment. Indeed, as seen in 
Fig. |l^, displaying the quasi-particle routhians characteristic of a well-deformed prolate g- 
band in ^''^Zr, the proton lgg/2 crossing frequency is higher than the neutron l/iii/2 crossing 
frequency. In most cases, the aligned (i//iii/2)^ bands are triaxial, with slightly negative 7 
values. 

For some Sr and Zr isotopes, at large spins (/~30-50) there appear prolate well- 
deformed minima with /32~0.32-0.38. The associated configurations are {vhii/2Y{T^g9/2Y 
and (z//iii/2)'^(z^«i3/2)^(7r5'9/2)^ (see Fig. [l^). 

The experimental and calculated angular momentum alignment for the well-deformed 
near-prolate yrast bands of ^'^''Zr and ^'^^Zr are presented in Fig. IT8|. Although our cal- 



culations are adiabatic and hence unable to reproduce the alignment pattern in the band 
crossing region, it is clearly seen that the first band crossing is associated with the i//iii/2 
alignment. 



2. Triaxial or near-oblate collective minima with -75°^j^—15° 

In the heavy Sr isotopes there appear excited minima at oblate shapes which give rise to 
rotational bands (g'-bands) at 7=-60°. Except for the A^=58 isotope (/?2=0.17), associated 
f32 deformations are typically ;^0.25 and increase slightly with increasing I. The oblate bands 
are predicted to lie at E*^l MeV above the yrast configuration. At higher spins, 1^20, 
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triaxial bands built upon several aligned quasi-particles [(i^/iii/2)^(7r5'9/2)^, {^hii/2)^{ngQ/2)'^, 
(i//iii/2)'^(i^ii3/2)"^(vr(yf9/2)^] emerge. Within these bands, P2 values vary smoothly, from /?2~0.3 
to ?s0.22. The triaxial many-quasi-particle bands extend over 1^24-46 for A^=60 and 62, 
/Ril8-52 for iV=64, and 7^30-52 for iV=66. They are highly excited (E*>2MeV above 
yrast). 

The heavy Zr isotopes with A^=58-62 are predicted to have oblate collective bands similar 
to those calculated for the Sr isotopes. For N=Q4 and 66, triaxial minima with negative 
7 deformations are predicted. The corresponding configurations are (i^/iii/2)^(7i"5'9/2)^ and 
{i^hn/2)'^{TTgg/2y, extending from 1^20 to 1^42. In ^°^Zr, a 8-q.p. (i^/iii/2)'^(i^«i3/2)H7r5'9/2)^ 
band is predicted at high spins, /~50. However, due to the large quasi-particle content, its 
collectivity is weak (/52~0.15, 7~-80°). The excitation energy of this band is 1-4 MeV above 
the superdeformed (SD) yrast band. 

The heavy Mo isotopes are predicted to have triaxial g-bands with negative 7 values. 
For A^=58 and 60, the triaxial minima are very shallow. The large excursions in the 7 
direction shown in Fig. |15| are due to this softness (or 7-instability) and hence should not 
be taken literally. The heavier Mo isotopes have better localized triaxial minima (/32~0.3, 
7~— 20°). Consequently, the deformation changes in their g-bands are small. The quasi- 
particle routhian diagram, characteristic of triaxial configurations in Mo and Ru isotopes, 
is presented in Fig. IS. Figure 20 shows the experimental and calculated angular momen- 
tum alignment for the yrast band of ^°^Ru. In contrast to the situation for prolate shape 
(Fig. |1^, here the 'n'gQ/2 alignment is expected to occur at a lower frequency. Because of 
the deformation softness, shape changes with quasi-particle alignment are more dramatic 
in the Mo isotopes than in their (more rigid) Sr and Zr isotones. In most cases, highly 
aligned bands have 7^—30°, and their collectivity gradually decreases with /. For exam- 
ple, the {i^hii/2)^{i'ii3/2Y{T!'g9/2)^ and (z^/iii/2)^('^^i3/2)"'^(7r5'9/2)'^ configurations calculated at 
1^40-50, have weakly deformed shapes with /52~0.12. 

In recent experimental work |2D|, transition quadrupole moments of high-spin states in 
the neutron-rich Mo isotopes have been measured using the Doppler-profile method. The 
observed reduction of quadrupole moments at J^IO has been interpreted in terms of a 
deformation change from axial shapes towards triaxial shapes with 7>0. Our calculations 
do not confirm this suggestion. According to the results shown in Fig. |15|, deformations 
of triaxial g-bands in, e.g., lO'^'ioe.iosj^Q ^^.^ rather stable, and it is only after the vhii/2 
alignment that some small reduction in (32 is predicted. In our opinion, the decrease in 
the transition quadrupole moments seen experimentally is due to the mixing between a g- 



band and a (2^/^11/2)^ band at 7:^10. Indeed, as noted in Ref. [^, the observed moments 



of inertia reveal a gradual increase for /ict;=0.2-0.4MeV, reflecting a large-interaction band 
crossing. There is a difference between the TRS calculations presented in Ref. and these 
in our study. Namely, the previous study predicts positive 7 deformation for the aligned 
{vhiii2f band in ^^^Mo. In our opinion, this discrepancy is due to the different spin-orbit 
parameterizations employed in these two papers. As discussed in Sec. 0, the spin-orbit 
potential used in our study was optimized to reproduce the experimental shape transition 
near iV=56. The change in the spin-orbit strength implies a different position of the 1/111/2 
neutron shell and, consequently, gives rise to a different predicted deformation pattern. 

The neutron-rich Ru isotopes have 7-soft ground states which are stabilized by rotation 
(see Fig. O). For iV=60-66, the g-bands correspond to triaxial shapes with negative values 
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of 7. (The positive-7 minima never are predicted to be yrast at low and medium spins.) 
The deformation changes in ^°^Ru due to and ahgnments are predicted to be 

rather small, with the g-band, (1^^11/2)^ band, and (z^/iii/2)^(7r5'9/2)^ band having similar 
deformations, (/?2, 7)~(0.25, —20°). The strongly triaxial {yhii/2)'^{j^g9/2Y configuration in 
^°^Ru becomes yrast at 1^40, but it is quickly crossed by a SD band. The nuclei ^°^'^^°Ru 
behave very similarly. Their g-bands correspond to deformation (/32, 7)~(0.28,— 20°). The 
vhiij2 alignment, the TTgg/2 alignment, and the second i^/iii/2 alignment produce triaxial 
shapes with /52~0.2, 7^—45°. The z//iii/2 crossing is calculated to occur at ^cc;~0.35MeV 
for both nuclei. This is slightly lower than the observed crossing frequency, ^cc;exp~0.4MeV 
(see Figs. |2g and 0). 

Properties of heavy Mo and Ru isotopes are often discussed in terms of triaxial degrees of 
freedom. Experimentally, the 2^ states lie at fairly low excitation energies (~700-800 keV), 
indicating the importance of 7 deformation. For ^°'^'^^°Ru, the B{E2) branching ratios mea- 
sured in the 7-band at low spin are consistent with that of a rigid triaxial rotor [|I^]. However, 
it has also been pointed out that at higher spins the analysis with the rotational-vibrational 
model gives a better agreement 0. More recent calculations based on the generalized col- 
lective model suggest that the experimental result is consistent with the assumption of 
a 7-soft potential with a triaxial minimum. It should also be pointed out that the CHFB 
calculations with the pairing-plus-quadrupole Hamiltonian |83,84] yield triaxial shapes, and 



that in the interacting boson model analysis of Ref. |76] these nuclei are calculated to be 



near the 0(6) limit rather than the SU(3) limit, thus indicating 7 softness. 



3. Triaxial minima with (32=0.20-0.35, 15°<'y<45° 

Triaxial minima with 7>0 are predicted at low angular momentum in the Mo and Ru 
isotopes. At 1=0 these minima are intrinsically equivalent to the triaxial ground-state 
configurations with 7<0. However, due to the large moment of inertia at 7<0, the latter 
bands are usually favored by rotation. In most cases, the 7>0 bands have {f32,l)^ (0.3, 
15°), and the deformation decreases to {^2,1) ~(0.2, 30°) in the corresponding aligned 
configurations. An exception is the nucleus ^"^^Mo, whose yrast line has {(32, l) ~ (0.2, 30°) 
up to 1^25 (the negative-7 minimum quickly disappears). 

At intermediate angular momentum (/?«24-30, huj ^0.6-0. 7 MeV) in 96,98, i02gj, ^^"^Zr 
appear well-deformed minima appear with (/32,7)~ (0.2, 30°). These correspond to the 
aligned (z^/iii/2)'^(7r5'9/2)^ configuration and become yrast or near yrast at 1^24. 

At very high spins (7^50-70, huj^l.0-1.5 MeV) in 98,ioo,io4,io62r and io2,io6,io8j^q strongly 
deformed triaxial minima with /32=0.28-0.4 and 7~30° are predicted. These bands involve 



the 7r/iii/2 and i'ii3/2 intruders (see Figs. p!4| and p!5D . 



4. NoncoUective configurations 

Noncollective many-quasi-particle states (7=60°) correspond to the optimal (stretched) 
shell-model configurations with the single-particle angular momenta of all decoupled nu- 
cleons aligned along the axis of the total angular momentum. Often, the stretched states 
correspond to termination points of collective rotational bands. The many-quasi-particle 
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bands with small (32 deformations and/or 7- values close to the 7=60° limit (seen in Figs. 



16) are candidates for such terminating structures. (For recent calculations of terminating 



configurations in this mass region, see Ref. [p5| .) 

In all nuclei considered, the TRS calculations predict that the noncoUective configurations 
approach yrast in the spin range 1=30-50. The corresponding deformations increase with 
increasing /, approaching /52~0.3-0.35 at 1^60. Since the main focus of our paper is on 
collective rotation and shape effects, and since the adiabatic TRS method used does not allow 
for a detailed description of optimal states, we do not discuss noncoUective configurations 
quantitatively. 



5. Superdeformed near-prolate minima with [32^0.40 

The nuclei from the heavy-Zr region are excellent candidates for finding SD bands. In our 
calculations, SD minima have been predicted in ^°^~^^°Ru, ^''°~^°^Mo, io2-io62j^^ g^^^j ^°°Sr. 
These minima approach yrast at high rotational frequencies, typically at /icu^O.S MeV. Since 
at these high angular momenta pairing correlations are less important ||105|| , intrinsic config- 



urations of SD states are well characterized by the intruder orbitals carrying large principal 
oscillator numbers M (high-A^ orbitals) | p.07] , |115| - |117|| . In this mass region, these are the 
Af=5 proton and Af=6 neutron states. Because of their large intrinsic angular momenta 
and quadrupole moments, high- A/" orbitals strongly respond to the Coriolis interaction and 
to the deformed average field. Consequently, their occupation numbers are good indicators 
of rotation and deformation properties of SD bands. 

Figure ^ shows the single-particle routhian diagram representative of discussed SD con- 
figurations. The single-particle gaps in the routhian spectrum indicate the best candidates 
for SD bands: at low spins, the combination of the Z=42 and A^=58 gaps leads to ^°°Mo, 
while the Mo, Tc, and Ru isotopes with N approaching 68 are expected to favor SD high 
spin bands. This observation is fully corroborated by our calculations. Indeed, the z/6°7r5° 
and z/6^7r5^ SD bands in ^°°Mo and ^^°Ru, respectively, are predicted to appear near yrast 
already at 1^22. 

For the heavy Sr isotopes, SD minima are calculated only for ^°°Sr, and only at very high 
spins, 1^64-70. The SD band contains two M=Q neutrons and one J\f=5 proton; hence the 
corresponding configuration is labeled z/6^7r5^. This configuration is also predicted in ^°^Zr 
and ^°^Zr. The corresponding SD bands behave rather smoothly with uj up to the highest 
spins considered. In contrast, for ^°^Zr the SD yrast line is less regular. Here, the i/6^7r5° 
band is crossed at ^a;i^l.l MeV by the more strongly deformed z/6^7r5^ band. At still higher 
spins, there is a transition to an even more deformed structure, z/6^7r5^. The corresponding 
kinematic moments of inertia, J^^^\ extracted from the TRS calculations are displayed in 



Fig. 22. It is seen that the moments of inertia in all z/6^7r5^ bands decrease gradually with 
frequency. The increased J'^^^ in ^°^Zr is due to the two band crossings discussed above. 

At low spin, the SD minimum of the A=58 isotones of Mo and Ru (^°°Mo and ^°^Ru) 
can be associated with the configuration i/6°7r5°. At hu^l MeV, the configuration z/6^7r5^ 
becomes lower in energy (see Fig. |22|). It is noted that ^°°Mo is among the few nuclei, 
together with ^°^Mo and ^^°Ru, in which the SD minimum appears at particularly low spins, 
i.e., at 1^20 (the corresponding SD band is calculated to appear ?^i3MeV above yrast). As 
discussed above, the microscopic reason for this structure can be explained in terms of the 
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SD shell gaps at Z=42 and N=58. The nuclei 102,104,106^^ g^^^^ i06,i08j^^ predicted to 
have rather regular SD bands associated with the z/6^7r5^ and i/6^7r5^ configurations. On the 
other hand, the SD yrast line of ^''^Mo and ^'''^Ru is expected to be less regular, with several 
consecutive band crossings. The best candidates for observing high-spin SD bands in this 
mass region are ^^°'^^^Ru. 



V. CONCLUSIONS 

The macroscopic-microscopic method based on the deformed Woods-Saxon potential was 
applied to neutron-rich nuclei from the A^lOO mass region. Thanks to a number of experi- 
mental developments, one is now able to obtain spectroscopic information on these systems, 
and more data should become be available in the near future, due to the large gamma-ray 
arrays such as GammaSphere or EuroBall, and the advent of neutron-rich radioactive beams. 

The main focus of this study is on the deformation properties and shape coexistence 
effects at low and high angular momentum. High-spin calculations were carried out for 
the vacuum configurations (i.e., positive-parity, even-spin) in even-even Sr, Zr, Mo, and Ru 
isotopes with 58<A^<66. 

The single-particle model used was previously optimized to reproduce the transition 
from spherical to deformed shapes at iVRi58. The calculated ground-state deformations for 
heavier systems are in general agreement with experimental data. 

The infiuence of particle number projection has been investigated. In general, the po- 
tential energy surfaces calculated with the PNP method have better developed deformed 
minima. A replacement of the energies at spin zero, E'^^^ in Eq. (|^), calculated within the 
BCS by those calculated with PNP, has only minor consequences for predicted high-spin 
properties. 

At low and medium spins, the alignment pattern is governed by an alignment of the i'hu/2 
and 7r(yf9/2 quasi-particles. The associated deformation changes are explained in terms of the 
shape polarization exerted by aligned high-j quasi-particles. According to the calculations, 
the first band crossing in the neutron-rich nuclei considered can be associated with the 
alignment of the l/;.ii/2 neutron pair. The frequency of the second bandcrossing, due to 
the 15^9/2 proton alignment, is rather high in the well-deformed Sr and Zr isotopes, but it 
becomes comparable with the neutron crossing frequency in the triaxial Mo and Ru isotopes. 

At higher angular momenta, the high-A/" intruder orbitals, namely A/'=5 protons and 
A/'=6 neutrons, become important. They are occupied in superdeformed structures and in 
some highly-triaxial configurations. 

From the point of view of nuclear quadrupole collectivity, probably the most interesting 
nuclei in the neutron- rich A^^lOO mass region are systems near ^"^Mo and ^°^Ru, which are 
predicted to have stable triaxial shapes (7~-30°) in 0-q.p., 2-q.p., and 4-q.p. configurations. 
Although the question of whether they are 7-soft or 7-deformed at low spins has not yet 
been settled, these systems seem to form a very good testing ground for theoretical models of 
nuclear triaxiality. In particular at higher spins, where the triaxial minima are predicted to 
be deeper, the shape with 7f»-30° can give rise to interesting selection rules associated with 
the effective C4 symmetry of the Hamiltonian | |118| ]. Finally, the most favorable candidates 



for super deformation in this mass region are i°°Mo and ^^^'^^^j^u in which the SD minima 
are predicted at particularly low spins. 
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TABLE L Calculated equilibrium deformations ^2,^4, and of the ground-state and excited 
minima in neutron-rich even-even Se-Cd isotopes. The energy differences between the local minima, 
and the corresponding quadrupole moments are also shown, together with experimental quadrupole 
moments |109|. 
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3.62 
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0.13 
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4.86 


3.68 
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0.33 
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-2.40 


A 1 f 

4.15 




70 


-0.25 
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4.02 
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-0.01 
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-1.42 
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3.49 
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0.62 
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2.91 
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-0.03 


-0.01 


0.25 


0.04 


-0.01 


0.45 


-2.40 


3.40 




64 


-0.24 


-0.05 


0.00 


0.26 


0.03 


-0.01 


0.15 


-2.45 


3.54 


3.22 


66 


-0.24 


-0.07 


0.00 


0.26 


0.01 


-0.01 


-0.36 


-2.47 


3.41 


3.34 


68 


-0.25 


-0.09 


0.00 


0.25 


-0.02 


0.00 


-0.88 


-2.51 


3.30 


3.36 


70 


-0.25 


-0.09 


0.00 


0.24 


-0.03 


0.00 


-1.19 


-2.52 


3.17 




72 


-0.24 


-0.10 


-0.01 


0.22 


-0.04 


0.00 


-1.43 


-2.41 


2.89 




74 


-0.21 


-0.11 


-0.03 


0.17 


-0.02 


0.00 


-1.58 


-2.12 


2.21 




76 


-0.18 


-0.11 


-0.04 


0.08 


0.00 


0.00 


-1.31 


-1.90 


1.11 












0.51 


0.04 


0.05 


5.08 




8.05 





TABLET. (Continued) 
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MeV 
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0.21 


2.15 


58 








0.14 


0.03 


-0.01 






1.89 


2.32 


60 








0.18 


0.03 


-0.01 






2.42 


2.57 


62 


-0.23 


-0.01 


-0.01 


0.21 


0.04 


-0.01 


1.00 


-2.60 


2.97 


2.76 


64 


-0.24 


-0.03 


-0.01 


0.22 


0.03 


-0.01 


0.74 


-2.63 


3.11 


2.96 


66 


-0.24 


-0.05 


0.00 


0.21 


0.00 


0.00 


0.37 


-2.64 


2.91 


2.52 


68 


-0.24 


-0.07 


0.00 


0.21 


-0.02 


0.00 


0.00 


-2.66 


2.82 


1.84 


70 


-0.24 


-0.08 


0.00 


0.20 


-0.03 


0.00 


-0.20 


-2.65 


2.77 


2.40 


72 


-0.23 


-0.08 


-0.01 


0.18 


-0.03 


0.00 


-0.24 


-2.54 


2.47 




74 


-0.19 


-0.09 


-0.03 


0.14 


-0.02 


0.00 


-0.16 


-2.12 


1.91 




76 


-0.13 


-0.07 


-0.03 


0.12 


-0.01 


0.00 


-0.07 


-1.56 


1.61 












0.51 


0.03 


0.06 


7.52 




8.40 




78 


0.00 


0.00 


0.00 


0.54 


0.03 


0.05 


-9.25 


0.03 


9.07 




48 56 








0.00 


0.00 


0.00 






0.00 




58 








0.03 


0.00 


0.00 






0.43 


2.03 


60 


-0.04 


0.00 


0.00 


0.12 


0.01 


0.00 


0.24 


-0.50 


1.66 


2.08 


62 


-0.07 


0.01 


0.00 


0.15 


0.00 


0.00 


0.42 


-0.90 


2.07 


2.13 


64 


-0.09 


0.01 


0.00 


0.17 


0.01 


0.00 


0.52 


-1.18 


2.49 


2.26 


66 


-0.10 


0.00 


0.00 


0.18 


-0.01 


0.01 


0.57 


-1.34 


2.61 


2.35 


68 


-0.11 


-0.01 


0.00 


0.17 


-0.03 


0.01 


0.56 


-1.37 


2.41 


2.37 


70 


-0.10 


-0.01 


0.00 


0.16 


-0.03 


0.01 


0.48 


-1.25 


2.31 




72 


-0.08 


-0.02 


0.00 


0.14 


-0.03 


0.01 


0.24 


-1.09 


1.96 




74 


-0.07 


-0.02 


-0.01 


0.09 


-0.02 


0.00 


0.07 


-0.90 


1.28 




76 








0.00 


0.00 


0.00 






-0.03 





TABLE I. (Continued) 



23 



PNP BCS 



z 


N 




/34 


7 


Qcal 


/?2 


/94 


7 


Qcal 


42 


58 


0.210 


0.011 


-32.726 


3.008 


0.162 


0.001 


-51.801 


2.103 




60 


0.300 


0.035 


0.000 


3.973 


0.243 


0.019 


-20.538 


3.547 




62 


0.316 


0.017 


-19.046 


4.737 


0.291 


0.021 


-14.055 


4.261 




64 


0.335 


0.016 


-19.695 


5.116 


0.308 


0.015 


-16.866 


4.619 




66 


0.328 


0.001 


-21.119 


5.008 


0.306 


0.005 


-18.690 


4.629 


44 


58 


0.231 


0.015 


-25.030 


3.541 


0.125 


0.013 


0.000 


1.606 




60 


0.252 


0.019 


-23.157 


3.940 


0.184 


0.020 


0.000 


2.460 




62 


0.275 


0.003 


-19.507 


4.271 


0.245 


0.003 


-18.232 


3.756 




64 


0.284 


0.014 


-19.875 


4.533 


0.263 


0.013 


-21.254 


4.185 




66 


0.283 


0.004 


-21.870 


4.550 


0.272 


0.012 


-23.660 


4.413 



TABLE II. Equilibrium deformations P2, 7, and (3^ and corresponding quadrupole moments 
(in eb) for the Mo and Ru isotopes calculated in the PNP (left) and BCS (right) variants. 
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FIGURES 



FIG. 1. Single particle neutron (top) and proton (bottom) levels of the deformed Woods-Saxon 
potential as a function of quadrupole deformation /?2. Positive (negative) parity states are indicated 
by solid (dotted) lines, and the spherical and deformed shell and subshell closures are indicated. 
Calculations are performed with the Woods-Saxon set of parameters of Ref. corresponding to 
the central nucleus ^'''^Zr. 

FIG. 2. Potential-energy curves for the even-even neutron-rich isotopes of Kr, Sr, and Zr ob- 
tained in the FRLD/PNP (left) and FRLD/BCS (right) variant of calculations. At each value of 
/32, the energy has been minimized with respect to deformations /?4 and (3%. 

FIG. 3. Same as in Fig. ^ except for the even-even neutron-rich isotopes of Mo, Ru, and Pd. 

FIG. 4. Total energy surfaces in the (/32, 7)-plane for 100.1022^^ LD/PNP (left) and 

LD/BCS (right) variant of calculations. At each mesh point (/32, 7) the total energy has been 
minimized with respect to fi^. The distance between thick contour lines is 1 MeV, while between 
the thin contour line it is 250 keV. 

FIG. 5. Same as in Fig. | except for I04,i06 2r. 

FIG. 6. Same as in Fig. | except for ioo,io2j^q^ 

FIG. 7. Same as in Fig. | except for lO^'^o^Mo. 

FIG. 8. Same as in Fig. | except for i^^'^o^Ru. 

FIG. 9. Same as in Fig. | except for i°^'^i°Ru. 

FIG. 10. Total routhian surfaces in the (/32, 7)-plane for the (7r=-|-, r=l) quasi-particle vacuum 
configuration of ^'^'^Sr at four values of rotational frequency: huj=Q.2i^ 0.6, 0.9, and 1.2 MeV. At 
each (/?2, 7) point the total routhian has been minimized with respect to hexadecapole deformation 
/34. The distance between thick contour lines is 1 MeV, while between the thin contour line it is 
250 keV. The angular momentum values at local minima are indicated. 



FIG. 11. Same as in Fig. |10| except for ^'^^Zr. 



FIG. 12. Same as in Fig. 110 except for ^'^^Ru 



25 



FIG. 13. Summary of calculated equilibrium deformations of yrast and near-yrast (7r=+, r=l) 
bands in 96,98,ioo,i02,i04gj,^ The rotational frequency is varied from huj=0 to ^u;=1.5MeV, with 
steps of Aa;=0.07 MeV/?i. If the two minima corresponding to the TRS's calculated at lv and 
co+Auj have the same intrinsic configuration, and the change in equilibrium deformation is small, 
A(/32 cos 7) < 0.1, they are assumed to belong to the same rotational band and they are connected 
by a solid line. Whenever feasible, the direction of increasing uj is marked by an arrow. The symbols 
indicating the intrinsic configuration are put at the lowest and highest frequency at which the band 
is calculated, and the corresponding spin range is shown by numbers. The legend displays band 
labels. Superdeformed bands and high-spin bands with 7>0 are classified according to the number 
of occupied high-AA intruder levels {J\f = 6 and 5 for neutrons and protons, respectively). Other 
band structures are classified according to the number of aligned high-j quasi-particles: i^/iii/2) 
TTfifg/a, and ^1^/2- 



FIG. 14. Same as in Fig. O except for 98,ioo,i02,i04,i062j.. 



FIG. 15. Same as in Fig. Ill except for ioo,i02,i04,i06,i08]y[Q_ 



FIG. 16. Same as in Fig. |T| except for I02,i04,i06,i08,ii0j^^_ 

FIG. 17. Single quasi-particle routhian diagram for N=62 (top) and Z=40 (bottom) at /52=0.3 
and 7=0°. This diagram is typical for the well-deformed prolate bands in the heavy-Zr region. The 
parity and signature of individual levels are indicated in the following way: 7r=+, r=—i - solid 
line; 7r=+, r=+i - dotted line; 7r=-, r=—i - dot-dashed line; 7r=-, r=+i - dashed line. The thin 
line indicates the Fermi energy. 

FIG. 18. Experimental and calculated angular momentum alignment for yrast bands of 
well-deformed ioo,i022j._ Experimental data are taken from Ref. |19]. The diabatic bands have 



not been extracted in the calculation. Therefore, the band crossing region is not reproduced theo- 
retically. 

FIG. 19. Single quasi-particle routhian diagram for A^=64 (top) and Z=44 (bottom) at 
/52=0.23 and 7=-30°. This diagram is typical for the triaxial bands in the heavy-Mo region. The 



line characteristics have the same interpretation as in Fig. 17 



FIG. 20. Experimental and calculated angular momentum alignment for triaxial yrast line in 
^'^^Ru. Experimental data are taken from Ref. The diabatic bands have not been extracted in 
the calculation. Therefore, the band crossing region is not reproduced theoretically. 

FIG. 21. Single-particle routhian diagram for neutrons (top) and protons (bottom) at /?2=0.45 
and 7=0°, characteristic of SD configurations in the heavy Zr region. The line convention is the 



same as in Fig. 17 
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FIG. 22. Calculated kinematic moments of inertia versus rotational frequency for superde- 
formed bands in even-even nuclei around ^''^Mo. The symbols indicate the same configurations as 



in Figs. 13-16 
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